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Many inflating spacetimes are likely to violate the weak energy condition, a key assumption of 
singularity theorems. Here we offer a simple kinematical argument, requiring no energy condition, 
that a cosmological model which is inflating - or just expanding sufficiently fast - must be incomplete 
in null and timelike past directions. Specifically, we obtain a bound on the integral of the Hubble 
parameter over a past-directed timelike or null geodesic. Thus inflationary models require physics 
other than inflation to describe the past boundary of the inflating region of spacetime. 
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/. Introduction. Inflationary cosmological models jJSQ 
are generically eternal to the future |4|, |£g. In these mod- 
els, the Universe consists of post-inflationary, thermal- 
ized regions coexisting with still-inflating ones. In co- 
moving coordinates the thermalized regions grow in time 
and are joined by new thermalized regions, so the comov- 
ing volume of the inflating regions vanishes as t — » oo. 
Nonetheless, the inflating regions expand so fast that 
their physical volume grows exponentially with time. As 
a result, there is never a time when the Universe is com- 
pletely thermalized. In such spacetimes, it is natural 
to ask if the Universe could also be past-eternal. If it 
could, eternal inflation would provide a viable model of 
the Universe with no initial singularity. The Universe 
would never come into existence. It would simply exist. 

This possibility was discussed in the early days of in- 
flation, but it was soon realized that the idea could 
not be implemented in the simplest model in which the 
inflating universe is described by an exact de Sitter space. 
More general theorems showing that inflationary space- 
times are geodesically incomplete to the past were then 
proved ||. One of the key assumptions made in these 
theorems is that the energy-momentum tensor obeys the 
weak energy condition. Although this condition is satis- 
fied by all known forms of classical matter, subsequent 
work has shown that it is likely to be violated by quan- 
tum effects in inflationary models 0, Such viola- 
tions must occur whenever quantum fluctuations result 
in an increase of the Hubble parameter H — i.e., when 
dH/dt > — provided that the spacetime and the fluc- 
tuation can be approximated as locally flat. Such upward 
fluctuations in H are essential for the future-eternal na- 
ture of chaotic inflation. Thus, the weak energy condi- 
tion is generally violated in an eternally inflating uni- 
verse. These violations appear to open the door again to 
the possibility that inflation, by itself, can eliminate the 
need for an initial singularity. Here we argue that this is 
not the case. In fact, we show that the general situation 



is very similar to that in de Sitter space. 

The intuitive reason why de Sitter inflation cannot be 
past-eternal is that, in the full de Sitter space, expo- 
nential expansion is preceded by exponential contraction. 
Such a contracting phase is not part of standard inflation- 
ary models, and does not appear to be consistent with 
the physics of inflation. If thermalized regions were able 
to form all the way to past infinity in the contracting 
spacetime, the whole universe would have been thermal- 
ized before inflationary expansion could begin. In our 
analysis we will exclude the possibility of such a con- 
tracting phase by considering spacetimes for which the 
past region obeys an averaged expansion condition, by 
which we mean that the average expansion rate in the 
past is greater than zero: 

ffav > 0. (1) 

With a suitable definition of H and the region over which 
the average is to be taken, we will show that the averaged 
expansion condition implies past-incompleteness. 

It is important to realize that the terms expansion and 
contraction refer to the behavior of congruences of time- 
like geodesies (the potential trajectories of test particles). 
It is meaningless to say that a spacetime is expanding at 
a single point, since in the vicinity of any point one can 
always construct congruences that expand or contract at 
any desired rate. We will see, however, that nontrivial 
consequences can result if we assume the existence of a 
single congruence with a positive average expansion rate 
throughout some specified region. 

While the past of an inflationary model is a matter of 
speculation, the attractor nature of the inflationary equa- 
tions implies that many properties of the future can be 
deduced unambiguously. According to the standard pic- 
ture of inflation, all physical quantities are slowly vary- 
ing on the scale of if -1 . In the vicinity of any point P 
in the inflating region, we can choose an approximately 
homogenous, isotropic and flat spacelike surface which 
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can serve as the starting point for the standard analy- 
sis of stochastic evolution . A simple pattern of ex- 
pansion is established, in which the comoving geodesies 
x = const in the synchronous gauge form a congruence 
with 77 > -J (87r/3)Gpoj where po is the minimum energy 
density of the inflationary part of the potential energy 
function. This congruence covers the future light cone 
of P. While large fluctuations can drive H to negative 
values, such fluctuations are extremely rare. Once in- 
flation ends in any given region, however, many of the 
geodesies are likely to develop caustics as the matter 
clumps to form galaxies and black holes. If we try to 
describe inflation that is eternal into the past, it would 
seem reasonable to assume that the past of P is like the 
inflating region to the future, which would mean that a 
congruence that is expanding everywhere, except for rare 
fluctuations, can be defined throughout that past. 

For the proof of our theorem, however, we find that it is 
sufficient to adopt a much weaker assumption, requiring 
only that a congruence with H av > can be continu- 
ously defined along some past-directed timelike or null 
geodesic. 

In Section II, we illustrate our result by showing how 
it arises in the case of a homogeneous, isotropic, and spa- 
tially flat universe. In the course of the argument we shed 
some light on the meaning of an incomplete null geodesic 
by relating the affine parameterization to the cosmologi- 
cal redshift. In Section III we present our main, model- 
independent argument. In Section IV we offer some re- 
marks on the interpretation and possible extensions of 
our result. 

77. A simple model. Consider a model in which the metric 
takes the form 



ds 2 = dt 2 - a 2 (t)dx : 



(2) 



We will first examine the behavior of null geodesies, and 
then timelike ones. 

From the geodesic equation one finds that a null 
geodesic in the metric ©, with affine parameter A, obeys 
the relation 



dX oc a(t) dt. 



(3) 



Alternatively, we can understand this equation by consid- 
ering a physical wave propagating along the null geodesic. 
In the short wavelength limit the wave vector is tan- 
gential to the geodesic, and is related to the affine pa- 
rameterization of the geodesic by fc' 1 oc dx^ /d\. This al- 
lows us to write dX oc dt/uj, where w = k is the physical 
frequency as measured by a comoving observer. In an ex- 
panding model the frequency is red-shifted as lj oc l/a(t), 
so we recover the result of Eq. (0 . 

From Eq. ©, one sees that if a(t) decreases suffi- 
ciently quickly in the past direction, then J a(t) dt can be 
bounded and the maximum affine length must be finite. 



To relate this possibility to the behavior of the Hubble 
parameter 77, we first normalize the affine parameter by 
choosing dX — [a(t)/a(tf)] dt, so dX/dt = 1 when t = tf, 
where tf is some chosen reference time. Using 77 = a/ a, 
where a dot denotes a derivative with respect to i, we can 
multiply Eq. by H(X) and then integrate from some 
initial time ti to the reference time tf. 



H(X)dX= -^T<1. 

X(ii) Ja(U) a {tf) 



(4) 



where equality holds if a(ij) = 0. Defining 77 av to be an 
average over the affine parameter, 
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\{t f )-X(ti) 



77(A) dX < 
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X(t f )-X(U) ' 



(5) 

we see that any backward-going null geodesic with 77 av > 
must have a finite affine length, i.e., is past-incomplete. 

A similar argument can be made for timelike geodesies, 
parameterized by the proper time t. For a particle of 
mass to, the four-momentum P M = mdx^/dr, so we can 
write dr = (m/E)dt, where E = P° is the energy of 
the particle as measured by a comoving observer. If we 
define the magnitude of the three-momentum p by p 2 = 
—gij P l P° , where i and j are summed 1 to 3, then E = 
For a comoving trajectory we have P % 
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and therefore dr = dt. For all others, p oc l/a(t) [17|, so 
we can write p{t) = [a(tf)/a(t)]pf, where pf denotes the 
value of the three- momentum p at the reference time tf. 
Combining all this, we find 



ftf fa(tf) 

/ H(r)dr= / 



to da 
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(6) 



where the inequality becomes an equality when a(i;) = 
Here Ef = 



V reV wncrc 



Vrcl 



to 2 and 7 = 1/ 
Pf/Ef is the speed of the geodesic relative to 



the comoving observers at time tf. Since the integral 
is bounded, the argument used for null geodesies can be 
repeated, with the average taken over proper time. 

777. The main argument. In this section we show that 
the inequalities of Eqs. (@J and © can be established in 
arbitrary cosmological models, making no assumptions 
about homogeneity, isotropy, or energy conditions. To 
achieve this generality, we need a definition of the Hubble 
parameter 77 that applies to arbitrary models, and which 
reduces to the standard one (H = a/ a) in simple models. 

Consider a timelike or null geodesic O ("the observer"). 
We assume that a congruence of timelike geodesies ( "co- 
moving test particles" ) has been defined along O 0] , and 
we will construct a definition for 77 that depends only on 
the relative motion of the observer and test particles. 
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FIG. 1: The observer's worldline O and two test particles. 



In order to motivate what we do, we first consider 
the case of nonrelativistic velocities in Minkowski space. 
Suppose that the observer measures the velocities of the 
test particles as a function of the time r on his own clock. 
At time t\ particle 1 passes with velocity u(t\), and at 
time T2 = t\ + At particle 2 passes with velocity u(t2). 
What expansion rate could he infer from these measure- 
ments? The separation vector between the positions f\ 
and r*2 of the two particles at ti is Ar = r \ — r*2 = u At, 
and its magnitude is Ar = |Ar |. Their relative velocity 
is Au = u\ — vii = — {du/dr)AT. The Hubble expansion 
rate is defined in terms of the rate of separation of these 
particles, which in turn depends on the radial compo- 
nent of their relative velocity, Au r = Au ■ Ar/Ar. The 
inferred Hubble parameter H is then 



H 



Ai 



Au - Ar 



u • [du/ dr) 



At 



(7) 



which will equal the standardly defined Hubble parame- 
ter for the case of a homogeneous, isotropic universe. The 
expression for H may be simplified to H = —d \nv re \/dT, 
where w ro i = \u\. The fact that H is the total derivative 
of a function of v re \ implies that the variation of w ro i is 
determined completely by the local value of H, even if 
the universe is inhomogeneous and anisotropic. 

We can now generalize this idea to the case of arbitrary 
velocities in curved spacetime. Let w M = dx^/dr be the 
four-velocity of the geodesic O, where we take r to be 
the proper time in the case of a timelike observer, or an 
affine parameter in the case of a null observer. Let u^(t) 
denote the four- velocity of the comoving test particle that 
passes the observer at time r. We define 7 = u„v v , so in 
the timelike case 7 may be viewed as the relative Lorentz 
factor — i^ cl ) between w M and lA In the null case, 

7 = dt/dr, where t is the time as measured by comoving 
observers, and t is the affine parameter of O. 

Consider observations made by O at times t\ and 
T2 = T\ + At, as shown in Fig. ^ where At is infinites- 
imal. Let Ar^ be a vector that joins the worldlines of 



the two test particles at equal times in their own rest 
frame. Such a vector is perpendicular to the worldlines, 
and can be constructed by projecting the vector —v^At 
to be perpendicular to u M : Ar M = — d''At + ju^At. 
In the rest frame of the observer this answer reduces 
for small velocities to its nonrelativistic counterpart, 
Ar M = (0,u' At). This is a spacelike vector of length 
Ar = |Ar M | = ^/j 2 — k At, where n = v^v^ is equal 
to 1 for the timelike case and for the null case. The 
separation velocity will be Aw M = —(Du^/dr) At, where 
D/dr is the covariant derivative along O. The covariant 
derivative allows us to compare via parallel transport vec- 
tors defined at two different points along O, and can be 
justified by considering the problem in the free- falling 
frame, for which the affine connection vanishes on O at 
t = T\. The radial component of this velocity will be 
Au r = — (Au^Ar^) /Ar, where the sign arises from the 
Lorentz metric. We define the Hubble parameter as 0] 



„ _ Aw r -v^Du^/dr) 
ti = — j — = t. 

Ar 7 — ^ 



(8) 



Since O is a geodesic, we have (Dv^/dr) = 0, and 
therefore 



7 — K dr v ' 



(9) 



where 



7" 



null observer (k = 0) 



Fh)= < l 7 +l (10) 

— In timclike observer [k = 1) 

2 7 — 1 

As in Section II, we now integrate H along O from some 
initial Tj to some chosen Tf : 

pHdr = F( lf )-F( li )<F( V ). (11) 

In the null case Ffaf) = jj 1 , which is equal to the value 
of dr/dt at tf, normalized in Section II to unity. 

Eq. 111(1 therefore reproduces exactly the results of 
Eqs. Q arid (j^J, but in a much more general context. 
Again we see that if -ff av > along any null or non- 
comoving timelike geodesic, then the geodesic is neces- 
sarily past-incomplete. 

IV. Discussion. Our argument shows that null and time- 
like geodesies are, in general, past-incomplete in infla- 
tionary models, whether or not energy conditions hold, 
provided only that the averaged expansion condition 
H av > holds along these past-directed geodesies. This 
is a stronger conclusion than the one arrived at in pre- 
vious work 8] in that we have shown under reasonable 
assumptions that almost all causal geodesies, when ex- 
tended to the past of an arbitrary point, reach the bound- 
ary of the inflating region of spacetime in a finite proper 
time (finite affine length, in the null case). 
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What can lie beyond this boundary? Several possibil- 
ities have been discussed, one being that the boundary 
of the inflating region corresponds to the be ginn ing of 
the Universe in a quantum nucleation event [l2j. The 
boundary is then a closed spacelike hypersurface which 
can be determined from the appropriate instanton. 

Whatever the possibilities for the boundary, it is clear 
that unless the averaged expansion condition can some- 
how be avoided for all past-directed geodesies, inflation 
alone is not sufficient to provide a complete description of 
the Universe, and some new physics is necessary in order 
to determine the correct conditions at the boundary [2(| . 
This is the chief result of our paper. The result depends 
on just one assumption: the Hubble parameter H has a 
positive value when averaged over the affine parameter 
of a past-directed null or noncomoving timelike geodesic. 

The class of cosmologies satisfying this assumption is 
not limited to inflating universes. Of particular interest is 
the recycling scenario |l4| , in which each comoving region 
goes through a succession of inflationary and thermalized 
epochs. Since this scenario requires a positive true vac- 
uum energy p v , the expansion rate will be bounded by 
H m - m = yJ&ttGpv/?) for locally flat or open equal-time 
slicings, and the conditions of our theorem may be satis- 
fied. One must look carefully, however, at the possibility 
of discontinuities where the inflationary and thermalized 
regions meet. This issue requires further analysis. 

Our argument can be straightforwardly extended to 
cosmology in higher dimensions. For example, in the 
model of Ref. |l5j brane worlds are created in collisions 
of bubbles nucleating in an inflating higher-dimensional 
bulk spacetime. Our analysis implies that the inflating 
bulk cannot be past-complete. 

We finally comment on the cyclic universe model [Hi ] 
in which a bulk of 4 spatial dimensions is sandwiched 
between two 3-dimensional branes. The effective (3 + f )- 
dimensional geometry describes a periodically expanding 
and recollapsing universe, with curvature singularities 
separating each cycle. The internal brane spacetimes, 
however, are nonsingular, and this is the basis for the 
claim [16} that the cyclic scenario does not require any 
initial conditions. We disagree with this claim. 

In some versions of the cyclic model the brane space- 
times are everywhere expanding, so our theorem imme- 
diately implies the existence of a past boundary at which 
boundary conditions must be imposed. In other versions, 
there are brief periods of contraction, but the net result of 
each cycle is an expansion. For null geodesies each cycle 
is identical to the others, except for the overall normal- 
ization of the affine parameter. Thus, as long as H av > 
for a null geodesic when averaged over one cycle, then 
H av > for any number of cycles, and our theorem would 
imply that the geodesic is incomplete. 
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